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1 Utilizing renormalons 

Perturbative QCD relies on factorization. By this one implies that an observ- 
able that depends on at least one hard scale Q can be expanded in powers 
and logarithms of A/Q, where A is the intrinsic QCD scale. At least up to 
some order in 1/Q, one must also be able to factor a 'short-distance' part from 
long-distance contributions, which are independent of the details of the hard 
process. At leading order in 1/Q, the long-distance contribution can be absent, 
like in e + e~ annihilation, or a product of parton distributions like for inclusive 
quantities in hadron-hadron collisions. Beyond leading order, little is known 
about power corrections, with exceptions like deep inelastic scattering. For 
event shapes or effects referred to as 'hadronization', it is not known how to 
express power corrections in terms of operators and asymptotic states. 

The renormalon approach to power corrections uses the fact that the lead- 
ing term in the power expansion already indicates the existence of power cor- 
rections, because the perturbative expansion of its short-distance coefficient 
diverges. This 'renormalon' divergence occurs, because certain higher-order 
diagrams (the simplest being 'bubble graphs') contain many powers of loga- 
rithms of a loop momentum, which make these diagrams sensitive to large dis- 
tances.]] Summing a divergent series requires a prescription. The prescription- 
dependence is best captured by the ambiguity of the Borel integral and takes 
the form (A/Q) a ln b Q/A, where infrared (IR) renormalons yield positive inte- 
gers for a. One can interpret the ambiguity as an ambiguity in separating long- 
and short-distances, much as the factorization scale dependence in separating 
coefficient functions and parton distributions. Since the physical observable is 
unambiguous, the ambiguity in defining perturbation theory must be matched 
by power corrections and this determines their Q-dependence, but not their 
magnitude, just as the evolution of parton distributions is perturbatively cal- 
culable, but not their initial values. 

The advantage of the method is that it is entirely perturbative, although to 
all orders. To some extent, the language inherited from studies of large-order 
perturbation theory is an unnecessary complication, since the set of diagrams 
that leads to a divergent series really only probes the IR sensitive regions of 
low-order skeleton-like graphs. It would be desirable to classify these regions 
systematically by extending standard methods of perturbative factorization El 
that identify logarithmic IR divergences to subleading, power-like IR sensitiv- 
ity. Meanwhile, most calculations are done in the formal Nf — > oo-limit, which 
selects diagrams with one chain of fermion loops at leading order. Provided 

b The logarithms also enhance the sensitivity to distances much shorter than 1/Q. The 
corresponding ultraviolet renormalons will not be discussed here. 
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the phase space of a cut fermion loop is integrated unweighted, the structure of 
power corrections inferred from renormalons in this approximation is equivalent 
to calculating the low-order diagrams with finite gluon mass A and Interpreting 
non-analytic terms in the small- A 2 expansion as power correctionsB 

The method has its limitations, precisely because it is purely perturbative. 
To go beyond classifying the expected power corrections for each observable 
separately, one needs additional assumptions, which do not follow from pertur- 
bative considerations alone, such as universality of power corrections, to relate 
different observables. And, of course, the IR sensitivity of Feynman diagrams 
might not exhaust all possibilities for power corrections. 

In the past two years, these ideas have been applied to observables that 
do not admit an operator product expansion, such as event shapes in e + e _ 
annihilation, jet observables and the Drell-Yan cross section. In these cases, IR 
renormalons provide genuinely new information about power corrections. This 
talk gives a somewhat qualitative overview and emphasizes the outstanding 
issues. 

2 Power corrections to event shape variables 

Event shapes are constructed from IR safe weights on hadronic final states 
in e + e~ collisions. Thrust, for example, is defined as T = max^(^ i \pi ■ 
n\)/ J"^ \pi\. The theoretical prediction is computed in terms of parton mo- 
menta, while hadron momenta are measured. Matching partons and hadrons 
is dealt with as a hadronization correction, which is obtained from Monte Carlo 
programs and accounted for in determinations of a s from event shapes. The 
fragmentation models built into Monte Carlo programs lead to hadronization 
corrections that vary like 1/Q with the cms energy Q. 

Theoretically one considers hadronization as a soft parton phenomenon 
that takes over from the parton shower at a certain typical hadronic scale 
fih. This separation scale is not uniquely fixed and the boundary between 
perturbation theory and hadronization is vague. Thus, probing the boundary 
of perturbation theory with renormalons may tell us more about hadronization. 
Event shape variables have been computed to thls^end both with a finite mass 
gluon in the lowest order gluon emission diagramHut] and in the approximation 
of a single chain of fermion loops til, in which case the region of small invariant 
mass of the qq pair in a cut fermion loop is the important one. The two methods 
are not equivalent in this case, because the invariant mass distribution depends 
on how each particular event shape weights the qq phase space. 

The calculation of the average (1 — T) with finite gluon mass leads to a 
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I/Q-correction 



(1-T) ]1/Q =K.-. 



(1) 



The l/Q-correction comes only from the two-jet region T — > 1, when the gluon 
momentum becomes very small. Multiple gluon emission diagrams modify the 
constant K, but only if all gluons are soft, so that again T — ► 1. Consequently, 
if the two-jet region is excluded from the average over T, we expect a smaller 
hadronization correction, 



The same conclusion applies to the heavy jet mass average (M%). The DEL- 
PHI collaboration has reanalyzed!] the energy dependence of event shapes by 
adding l/Q and l/Q 2 terms to the next-to-leading order perturbative expres- 
sion (evaluated at scale fi = Q) and by fitting the coefficients of the power 
corrections to the data. No hadronization correction from Monte Carlo pro- 
grams is applied. Some of the results are reproduced in Tab. [I] and agree qual- 
itatively with the above predictions. The energy-energy correlation (EEC) is 
predicted □ to have l/Q-corrections at all angles, because the soft gluon re- 
gion contributes at all angles. It is also easy to see that the three-jet rates 
i?3 computed from the JADE clustering algorithm have l/Q-corrections, while 
the Durham algorithm should have only l/<5 2 -corrections, because it weights 
the region of soft gluons quadratically with their energy rather than linearly. 
The DELPHI analysis does not have enough data points to test this prediction 
for the Durham algorithm. 

The examination of expected power corrections provides some guidance 
to selecting 'good' event shapes, the good ones being those less sensitive to 
hadronization. To go further, one has to make the stronger assumption that 
hadronization corrections in the two-jet limit are universaloCl This implies 
that although the constants K above are not calculable for any observable, 
their ratio for different observables is calculable, because multiple soft parton 
emission modifies K in a universal way. Thus, fitting a l/Q-correction to one 
observable would determine the hadronization parameter once and for ever. 
The assumption of universality could be justified diagrammatically if an event 
shape variable did not resolve the soft parton kinematics, which in fact it does. 
□ For example, in the two-jet region 1 - T w (M%/Q 2 ) + {Mf/Q 2 ), where Mi is 
the light jet mass. If we now consider the diagram where a single emitted gluon 
splits into a qq pair, we find 1 — T = Mj^/Q 2 if both quarks are emitted into 
the same hemisphere, and 1 — T = 2M 2 /Q 2 if they are emitted into opposite 



(1 ~ ^ 1 )|l/Q,0.5<T<0.8 
(! _ T )\l/Q,al\T 



cx a s (Q) . 



(2) 
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Table 1: Fits to the Q-dependence of event shape variables. a s (Mz), the coefficient of 
a 1/Q-term, C\, and l/Q 2 -term, C'2 (not quoted) are fitted to obtain the second entry for 
each observable. For the first entry C2 is fixed to zero. 



Observable 


Ci/GeV 


a 3 {M z ) 


<1-T) 


0.82 ±0.07 
0.83 ±0.20 


0.123 ±0.002 
0.122 ±0.004 


O^l-T) 


0.37 ± 0.05 
0.20 ± 0.05 


0.121 ±0.008 
0.134 ±0.003 


{Ml/E 2 vls ) 


0.54 ± 0.08 
0.75 ±0.26 


0.121 ±0.002 
0.116 ±0.006 


J^dM h {Ml/El is ) 


-0.01 ±0.03 


0.123 ±0.000 


{M 2 d /El is ) 


0.19 ±0.04 
0.10 ±0.05 


0.094 ± 0.003 
0.097 ± 0.003 


f° Q 5 5 d cos 9 EEC 


1.68 ±0.05 
0.27 ±0.23 


0.115 ±0.002 
0.137 ±0.004 


Ri (ycut = 0.08) 


0.44 ±0.15 
-3.59 ±0.55 


0.107 ±0.001 
0.123 ±0.002 


RsiVout = 0.04) 


-0.67 ±0.49 
-2.53 ±3.15 


0.126 ±0.004 
0.137 ±0.019 



hemispheres. Thus, there is no unique relation between 1 — T and M^, even 
if both quarks are soft. -Universality could still hold in an approximate sense, 
if, as advocated in Ref . u, the strong coupling approaches a finite and not too 
large value in the infrared. In this case, the diagram just discussed is higher 
order in the IR coupling. In this scenario, the 1/Q-correction to should be 
smaller than for M^, because it arises only at second order. The small power 
correction for the average Mj = M? — Mf does not support this picture, 
although the small fit value for a s (Mz) indicates that the corresponding C\ 
in Tab. |] might not be too reliable. 

Eventually, universality should be subjected to experimental tests. In this 
respect, it would be interesting to obtain the coefficients C\ in Tab. |l| with 
u s fixed to a unique value. As a matter of principle, the power corrections 
obtained by renormalon methods are synonymous with large perturbative cor- 
rections in higher orders. If large coefficients are a practical concern, the 
divergent piece of the series should be separated and discarded, so that it is 
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absorbed into the power correction, leaving an unambiguous perturbative se- 
ries. A procedure of this sort has been proposed inQ and applied with some 
success to (1 — T), the average C-parameter and at,. Another question of 
importance for testing universality is to what extent the power corrections in 
Tab. [l] effectively parameterize higher order corrections in perturbation theory, 
which would in principle.be calculable, leaving a rather small 'true' hadroniza- 
tion correction. In Ref.tj it was argued that higher order corrections, summed 
up to the point where the series diverges, can well mimic the shape of a 1/Q- 
correction. An equivalent effect is obtained, if one expresses the second order 
perturbative prediction in terms of a s (Q*) with Q* ~ 0.1Q. Such a low scale 
is not unnatural for event shape observables, since they are dominated by the 
soft-collinear region, where the scale is set by the transverse momentum of 
the emitted gluon rather than Q. From this point of view, the question of 
whether universality holds is less important than determining the higher-order 
perturbative corrections or correct scale for each event shape. 

In principle, the universality assumption could also be invoked to relate two 
non-IR safe event shapes to each other. This would circumvent the difficulty 
of having to extract subdominant power corrections to test universality. 



3 Drell-Yan production and Sudakov resummation 



Drell-Yan (DY) production, apart from its phenomenological significance, is 
theoretically interesting, because one can kinematically realize the situation of 
a process with two hard scales. In the following, we consider the partonic DY 
cross section oby (after collinear subtractions) in the region Q ^> Q(l — z) ;§> 
A, where z = Q 2 /s, Q 2 being the mass of the DY pair and s the partonic cms 
energy. The second scale Q(l — z) can be identified with the energy available 
to parton emission into the final state. Since Q(l — z) <C Q, these partons 
are referred to as soft, although they are not soft in terms of the QCD scale 
A. Taking moments in z (roughly, this replaces 1/(1 — z) by N), one obtains 
two powers of In TV for each power of a s = a s (Q) , so that the actual expansion 
parameter of the hard cross section is a s In 2 N. Thus, in higher orders, one 
has two sources of large corrections, Sudakov logarithms, related to the scale 
Q(\ — z) and renormalon factorials, related to the scale A. One may ask how 
this complication affects the arguments that lead to the identification of power 
corrections through renormalons. 
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This question has been addressed in Refs.midO Starting from 

l >r Q 2 (i-zf 



2C F f, z N - 1 f dk 2 , , 
\na DY (N) = —±dz-^—— / -^a s {k x ), (3) 

Q 2 (l-z) 

which resums all leading logs InN (a s \nN) k in the DIS scheme, one finds E 
that the integral has a renormalon ambiguity of order NA/Q from the region 
of large z. However, the corresponding n\ occurs in far subleading logarithms, 
beyond the accuracy to which (||) was derived. Keeping all subleading loga- 
rithms essentially implies that the hard cross section is evaluated exactly. In 
Ref.lij this has been done in the approximation of a single chain, interpreted 
as an approximation to the logarithm of the cross section. Using the equiva- 
lence of this approximation with taking an explicit IR cut-off, we choose, for 
illustration, a lower cut-off /i on the emitted gluon energy. Omitting all terms 
that can not give rise to a /i/Q-correction (which allows us to ignore virtual 
corrections and collinear subtractions), one obtains instead of the right hand 
side of (|) 

l-2/*/Q Q 2 (l-z) 2 /4 

n J J k\ ^/(l- z) 2 -4k 2 ± /Q 2 



which reduces to the structure of (^) in the double-logarithmic, soft-collinear 
limit, if k± is set to zero in the square root. In this limit, consistent with the 
previous result, the integral contains a /i/Q-term in the expansion in small 
/I. However, the fcj_-term is crucial and can not be neglected precisely in the 
region z — * 1, where the fi/Q-teTin originates from. Keeping k± in the square 
root, one finds that l/Q-corrections are absent and that the leading power cor- 
rection is of order (N n/Q) 2 ■ The cancelation of the l/Q-term emphasizes that 
leading power corrections stem from soft gluons, but small angle (collinear) 
and large angle (kj_ ~ kg ~ Q(l — z)) emission are both important, til This 
might appear surprising, because large angle soft emission is usually considered 
suppressed, leading to angular ordering in parton cascades. But only logarith- 
mic enhancements of matrix elements cancel at large angles and no conclusion 
follows for power corrections. 

Returning to resummation of Sudakov logarithms, we conclude that there 
is no direct connection between resummation of In N terms and power correc- 
tions, which are 'buried' among an infinite number of subleading logs. Roughly 
speaking, this is so, because Q 3> Q(l — z) 3> A and Sudakov resummation 
is concerned only with the first inequality, while power corrections are asso- 
ciated with the smallest scale A. Nevertheless, one would like to formulate 
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the resummation procedure in such a way that it does not introduce stronger 
power corrections than required. Given the hierarchy of scales, it is useful to 
think about the problem from the effective field theory point of view, which 
deals with scales sequentially. Thus, one would first 'integrate out' momenta 
larger than Q(l — z) and deal with Sudakov resummation in this step. At that 
lower scale, we can consider power corrections, which then appear as the ratio 
A/(Q(l — z)) (NA/Q in moment space). 

The effective fields for fast quarks interacting with soft gluons are expressed 
as eikonal or Wilson lines. When N ^> 1, the hard DpJLYan cross section 
(omitting for simplicity collinear subtractions) factorizesliStj as &dy(N, Q) = 
Hdy(Q,^) Sdy{Q/N,h), where the scales Q and Q/N are separated. The 
'soft part' S satisfies a renormalization group equation in /i that can be used 
to sum logarithms in N, because S depends only on the single dimensionless 
ratio Q/(N/i). The solution to the RGE is expressed as 



ody = H DY (Q)- S DY (l,a s (Q/Nj) 




k 2 N 2 
r e ifc (ct s (fct))ln * 



+T DY (a s (kt)) 



(5) 



The three factors on the r.h.s. correspond to coefficient function, matrix el- 
ement and anomalous-dimension factor in the effective field theory language. 
The analysis of Ref. til shows that in the MS scheme the universal eikonal 
anomalous dimension and DY specific anomalous dimension Tdy are analytic 
functions for small a s , so that the exponential factor does not contain any 
renormalons (power corrections) at all. These enter only through the bound- 
ary conditions at the lower and higher scale and turn out to be (NA/Q) 2 as 
stated before. Notice also that when N > Q/A, the exponent becomes ill- 
defined. For such large moments, the language of power corrections looses its 
meaning (NA/Q ~ 1), the energy of soft gluons becomes of order A and any 
short-distance expansion fails. 

Whether the absence of 1/ Q-corrections for Drell-Yan production persists 
beyond the approximation of a single chain (one-gluon emission) , is an unsolved 
problem. In Ref. L3, the cancelation at leading order has been reproduced 
as a consequence of the KLN and Low theorems and it has been argued to 
hold at the level of two-gluon emission as well,— .On the other hand, in the 
language of Wilson lines, emphasized in Refs.ll3'Ej, a l/Q-correction could be 
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naturally accommodated by a certain operator constructed from Wilson lines, 
although it vanishes at leading order. If, as suspected in Ref. E3, Glauber 
gluons constitute a new potential source of l/Q-corrections at higher orders, 
the validity of the eikonal approximation and Wilson line treatment to power- 
like accuracy would have to be re-examined. In principle, the possibility exists 
liil that the cancelation of 1/Q terms occurs between (||) and terms dropped in 
(||), although there is no indication of it at leading order. Since, as explained 
above, the problem of resummation is disconnected from the problem of power 
corrections, one should be able to establish equivalence of the renormalon 
analysis for the DY cross section with the. analysis of power corrections in 
terms of multi-parton correlation functionstJ The present discussion indicates 
that the optimal language for problems with two scales might yet have to be 
found. 

4 The Sudakon: x-space vs. iV-space resummation 

As pointed out in Ref.0, factorial divergence, not related to renormalons, can 
appear when one converts resummed distributions in moment space (iV-space) 
back to 'real' space (x-space; note that x replaces z in the previous section). 
Typically, a physical quantity is given as an integral 



where a could be the partonic Drell-Yan cross section - in which case W is 
the parton luminosity and r = Q 2 / s - or the thrust distribution, or the lepton 
energy distribution in scmileptonic B — > X u lv decay, for example. Let x be 
a generic variable, such that x — ► 1 corresponds to the soft gluon region. If 
the weight function W constrains x to the region of large x, but such that 
1 — x is still large compared to A/Q, where Q is the generic hard scale, the 
perturbative expansion of er(r) contains the usual renormalon factorials and 
the corresponding poles in the Borel transform. 

Consider now the double logarithmic approximation with fixed coupling 
a s , in which all (a s In 2 N) k have been resummed in moment space. It is 
perfectly consistent to perform the inverse Mellin transformation to ir-space in 
the same approximation. Then a(x) « exp(ca s ln 2 (l — x)) with some constant 
c. Expanding cr(r) in a s , one finds the divergent series 



l 




(6) 
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~^F n (r)(4 C )"n!a» +1 (7) 

n 

independent of W, as long as it does not depend exponentially on x. In prac- 
tical cases c is such that this series diverges much faster than expected from 
renormalons. Continuing the tradition of misnomers, I call the corresponding 
singularity in the Borel plane a Sudakon pole. It is unphysical and appears, 
because in the process of resummation, one has dropped terms, which after 
integration over x give equally large contributions and cancel the singularity. 

On physical grounds, one expects Sudakov suppression, so c is negative. 
Then the integral in (Q) can just be done and there is no reason to re-expand 
it in a s . Even in this case, however, resummation has modified the analytic 
structure of the Borel transform, although the spurious pole is Borel summable. 
While c is indeed negative for the thrust distribution or lepton energy distri- 
bution in B decay, mentioned above, it is in fact positive for Drell-Yan produc- 
tion (and other hadronic collisions) in the conventional subtraction schemes, 
because the subtracted product of partonic structure functions shows stronger 
Sudakov suppression than the partonic Drell-Yan cross section. In this case, 
the integral in (Q) does not exist and must be defined by truncating the di- 
vergent expansion at its minimal term. If, as usual, one interprets the size 
of the minimal term as an uncertainty in defining the integral, one finds that 
this .uncertainty is of order (A/Q) l3o ^ 4c . As /?o/(4c) can be much smaller than 
one 113, this uncertainty is large. Alternatively, the integral can be defined 
by excluding the region of large x. Truncation of the series is equivalent to 
a cut-off in x that corresponds to excluding gluons with energy larger than 
A • (Q/A) 1-/3 °/ ( ' 2c ) from a perturbative treatment. This cut grows with Q and 
is much larger than the expected limit of order A for perturbative gluons. In 
fact, the energy cut corresponds to the position of the Sudakov peak of the 
x-distribution, if c were negative. 

Although the sketched procedure is consistent from the point of view of 
summing logarithms in 1—x, it is desirable to formulate resummation for inte- 
grated quantities as in ([|) such that factorial behaviour in .their expansions is 
consistent with the expected renormalon structure. In Ref.ED it is proposed to 
perform the inverse Mellin transform of the double-log resummed o'(N) back 
to x-space exactly, keeping all subleading logarithms. In this procedure, re- 
summation by itself does not introduce any factorial behaviour, as suggested 
also by the discussion of Sect. 3. Exact evaluation of the inverse Mellin trans- 
formation keeps exactly those subleading logarithms, which are necessary to 
cancel the unphysical Sudakon pole. This, together with the discussion of 
l/Q-corrections before, emphasizes that often it is not sufficient to perform 
resummations that are consistent to a certain logarithmic accuracy. Different 
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treatments of subleading logarithms can result in numerically important dif- 
ferent constant terms. Renormalon considerations can help to decide whether 
such terms are spurious or not. 
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